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Abstract 

We investigate the existence of bound states in a one-dimensional quantum system 
of N identical particles interacting with each other through an inverse square potential. 
This system is equivalent to the Calogero model without the confining term. The effec- 
tive Hamiltonian of this system in the radial direction admits a one-parameter family 
of self-adjoint extensions and the negative energy bound states occur when most general 
boundary conditions are considered. We find that these bound states exist only when 
N = 3, 4 and for certain values of the system parameters. The effective Hamiltonian for 
the system is related to the Virasoro algebra and the bound state wavefunctions exhibit 
a scaling behaviour in the limit of small inter-particle separation. 
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1. Introduction 



One dimensional many body quantum systems where the 2-body potential falls as the 
inverse square of the inter-particle distance exhibit rich mathematical structures and often lead 
to interesting physical applications [1-6]. This problem is equivalent to the Calogero model 
without the confining term. The spectrum of such a system is normally described in terms of 
scattering states alone ||. 

It is however known that 2-body systems on a half line containing an interaction term of 
the type where g is the 2-body coupling and x is the inter-particle distance, admit negative 
energy bound states for a certain range of parameter g J7|, These bound states owe their 
existence to self-adjoint extensions 11 of the corresponding Hamiltonian and have recently 



found physical applications |L0], |TT 



Another feature of these 2-body systems is that their corresponding Hamiltonians are re- 
lated to the Virasoro algebra ||12|| . Although the Hamiltonians themselves are not elements 



of the Virasoro algebra they do belong to the corresponding enveloping algebra. They are 
thus contained in the representation space of the Virasoro algebra. This fact has interesting 



implications for physics of black holes [11 



It is therefore natural to ask if these properties of the 2-body systems can be generalized to 
the iV-body case. In this Letter we provide a detailed analysis of this problem. We consider 
a system of iV identical particles of mass m interacting with each other through a potential of 
the type where x is the interparticle distance. The 2-body coupling g is assumed to be the 
same for all pairs of particles. The effective potential describing this system is found to be a 
function of g and N. The corresponding effective Hamiltonian is shown to admit self-adjoint 
extensions labelled by e tz where z G R. Different values of the parameter z leads to inequivalent 
quantizations of the system. The spectrum of the system is in general a function of g, N, and z. 
We show that for certain values of these parameters the system admits negative energy bound 
states. We also point out that the effective Hamiltonian for the iV-body system is contained 
in the enveloping algebra of the Virasoro algebra and thus have a well defined action on the 
representation space of the Virasoro algebra. The wave functions describing the negative energy 
bound states exhibit scaling behaviour in a certain regime. 

The organization of this paper is as follows. In Section 2, we study the iV-body system 
where the pairwise potential between any two particles is of the inverse square type and obtain 
the effective Hamiltonian for the system. In Section 3 we quantize the effective Hamiltonian and 
find the corresponding self-adjoint extensions. The negative energy bound states are obtained 
for certain values of the parameters of the problem. In Section 4 we indicate the connection of 
this system with Virasoro algebra and discuss the scaling behaviour of the bound states. We 
conclude in Section 5 with a summary and a outlook. 
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2. Effective Hamiltonian for the iV-body System 



Let us consider a system of N identical particles of mass m in one dimension whose Hamil- 
tonian is given by 

N q2 N i-1 



H E £3 2 + E E I 

i=l ux i i=2 j=l \ Xl 
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(2.1) 



where Xi is the coordinate of the i th particle. The interaction strength g between any two 
particles is assumed to be the same for all pairs of particles and units have been chosen such 
that 2mh~ 2 = 1. If g < — | the two-body system is unstable and the corresponding ground 
state energy is unbounded from below [13]. In what follows we shall assume that g > —\. 

We are interested in finding normalizable solutions of the eigenvalue problem 

Hip = Sip (2.2) 

especially when £ < 0. Following M, we consider the above eigenvalue equation in a sector of 
configuration space corresponding to a definite ordering of particles given by X\ > x 2 > ■ ■ ■ > 
xat. The translation- invariant eigenfunctions of the Hamiltonian H can be written as 



4> = y a+ H(r)P k (x), 



where 



N i-1 

nno* 

i=2j=2 
I N i-1 

iV i=2.7=l 

2 V 



Xj), 



2gY 



(2.3) 

(2.4) 

(2.5) 
(2.6) 



and Pk{x) is a homogeneous as well as translation- invariant polynomial of degree k which satisfy 
the equation 



N Q2 

E^ • - 



i=l 



dx 1 . 



N i-1 

EE 

i=2 j=l 



-{ 9 - 




\dxi 


dx jJ_ 



Pklx) 



0. 



(2.7) 



From Eqn. (2.3) we see that the condition a + | > must be satisfied for the wavefunction to 
be well behaved when Xi — > Xj. However, for the case a + | = 0, (2.1) evidently reduces to a 
free Hamiltonian of N particles. At present we shall consider only the case a + | > for which 
the wavefunction (2.3) vanishes at the limit Xi — > x 3 - 0. 

Substituting Eqn. (2.3) in Eqn. (2.2) and using Eqn. (2.7) we get 



dr 2 



r dr 



{21 
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where 



v = k+±(N-3) + ±N{N-l){a + ±). 



(2.9) 



We note that for A" > 3, v is a positive definite quantity. Let us define a function x(r) through 
the transformation 

l r , A 



(r ) = exp 



- J (1 + 2v)-dr x{r) = r ^ +u \(r). 



Substituting Eqn. (2.10) in Eqn. (2.8), we get 



dr'* 



+ 



v 



X(r) = Sx{r). 



(2.10) 



(2.11) 



As mentioned earlier, we are interested in solutions of Eqn. (2.2) where £ < 0. Writing £ = —E 
where E > 0, we can rewrite Eqn. (2.11) as 



Hx(r) = -Ex{r) 



where 



H 



d 2 



r-2 



(2.12) 



(2.13) 



v~ — j. H defines the effective Hamiltonian for the system and g is the effective 



dr 2 

and g = 
coupling. 

It may be noted that, the variable r in Eqn. (2.5) takes values in the positive real axis R + = 
[0, oo]. In fact, r represents the radial variable associated with — 1 number of translation- 
invariant 'Cartesian' coordinates (i.e. Jacobi coordinates), which are obtained from the Xj's 
through the elimination of center-of-mass coordinate H. Consequently, the operator H in Eqn. 
(2.13) defines the effective Hamiltonian in the radial direction for the A^-body system under 
consideration. By using Eqns. (2.3) and (2.10), the wavefunction of the A^-body system can 
be expressed in terms of the eigenfunction of the effective Hamiltonian as 



(2.14) 



We would next like to find the measure for which the eigenfunctions x{ r ) °f the effective 
Hamiltonian H would be square integrable. Let the angular variables associated with A^ — 1 
translation- invariant Cartesian coordinates be denoted by (i G [1, AT — 2]) 0. The relative 
distance between any two particles can be expressed in the form Xi — Xj = rXij(Q), where 
Xij(Q) is some function of these angular variables alone. Substituting this form of Xi — Xj 
into Eqn. (2.14) and also using the fact that Pk{x) is a translation- invariant homogeneous 
polynomial of degree k, if) can be factorized as ip — £,( r )F(£l) where 



f(r) = r^^+D+Mi+^r) (2.15) 

and F(Q) depends on the angular variables alone. Since if) is a normalizable eigenfunc- 
tion depending on ]V- 1 number of translation-invariant Cartesian coordinates, the integral 
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Jo°° £*( r )£( r ) r7V 2 dr must be finite and £(r) G L 2 [i? + ,r iV 2 dr]. Now, using Eqns.(2.15) and 
(2.9) we get, 

poo poo 

/ £*(r)t(r)r N - 2 dr = / x*(r)x(r)dr. (2-16) 
Jo Jo 

Thus / °° X*{ r )x{ r )dr is also finite and xi r ) £ £ 2 [-R + , <ir]. In the next Section we shall study 
the eigenvalue problem for the effective Hamiltonian. 

3. Self-adjoint Extension and Bound States 

The operator H of Eqn. (2.13) belongs to a general class of objects known as unbounded 
linear differential operators on a Hilbert space |§. We start by summarizing some properties 
of these operators which are needed for our purpose. 

Let T be an unbounded differential operator acting on a Hilbert space 7i and let (7, 5) 
denote the inner product of the elements 7, 5 G Ti. By the Hellinger-Toeplitz theorem ||, T 
has a well defined action only on a dense subset D(T) of the Hilbert space Ti. D(T) is known 
as the domain of the operator T. Let D(T*) be the set of <p G 7i for which there is an rj G Ji 
with (T£, 0) = (£, rj) V £ G D(T). For each such G D(T*) we define T*0 = 77. T* is called the 
adjoint of the operator T and D(T*) is the corresponding domain of the adjoint. 

The operator T is called symmetric or Hermitian if T C T*, i.e. if D(T) C D(T*) and 
T0 = T*0 V G D(T). Equivalently, T is symmetric iff (T<£, 77) = (0, Tr]) V 0, 77 G D(T). The 
operator T is called self-adjoint iff T = T* and D(T) = D(T*). 

We now state the criterion to determine if a symmetric operator T is self-adjoint. For this 
purpose let us define the deficiency subspaces K± = Ker(z =f T*) and the deficiency indices 
n±(T) = dim[i^±]. T is (essentially) self-adjoint iff (n+,n_) = (0,0). T has self-adjoint 
extensions iff n + = n_. There is a one-to-one correspondence between self-adjoint extensions 
of T and unitary maps from K + into K^. Finally if n + 7^ n_, then T has no self-adjoint 
extensions. 

We now return to the discussion of the effective Hamiltonian H. This is an unbounded 
differential operator defined in R + . H is a symmetric operator on the domain D(H) = {(f)(0) = 
0'(O) = 0, (f), (f)' absolutely continuous}. It is known that for g > |, H has deficiency indices 
(0,0) and is (essentially) self-adjoint on the domain D(H). For — | < g < |, H has deficiency 
indices (1,1) In the latter case is not self-adjoint on the domain D(H) but admits 

self-adjoint extensions. The deficiency subspaces K± in this case are one dimensional and are 
spanned by 

(f>+(r) = r^ijW(re^), (3.1) 
0_(r) = r5//( 2 )(re-^) (3.2) 

respectively, where -ff„'s are Hankel functions. The unitary maps from K + into K- are 
parametrized by e lz where z G R. The operator is self-adjoint in the domain D Z (H) which 
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contains all the elements of D(H) together with elements of the form + (r) + e* 2 0_(r). Different 
values of the parameter z lead to inequivalent quantizations of the operator H. 

As mentioned earlier, we are interested in the negative energy normalizable solutions of Eqn. 
(2.12). For g > | there are no such bound states ||. For each value of g such that — j < g < 3 
i.e. < v < 1, the solution of Eqn. (2.12) is given by 



4' 



X{r) 



B(VEr)2[J u (iVEr) - e lnu J_ V {W Er)\ 



(3.3) 



where J± u are Bessel functions and B is the normalization constant ||. In order to find the 
energy E, we use the fact that if H has to be self-adjoint, the eigenfunction x(r) must belong 
to the domain D Z (H). Note that in the limit r — > 0, 



r + e l 



smz/7r 



r(i + ^ 



— + " 



e 4 



r(i 



and 



B 



2e 2 



r i 



v) 



(3.4) 



(3.5) 



Thus if x( r ) ^ E) Z (H), then the coefficients of of r y+ 2 and r in Eqns. (3.4) and (3.5) must 
match. Comparing these coefficients we get 



E 



sin 



+ 3<) 



sm(f +7r|; 



(3.6) 



Thus we see that for every value of g within the range —j<g<j,H admits a single bound 
state with wavefunction given by Eqns. (3.3) and bound state energy given by £ = —E. 

We shall now show that the existence of the negative energy bound states for the original 
iV-body system depends on iV and the degree k of the polynomial Pk{x). To this end, let us 
recall that v must lie in the range < v < 1 for the negative energy bound states to exist. 
This together with Eqn. (2.9) gives 



o<fc + i(jv-3) + ^(jv-i)( a + i; 

Using Eqn. (3.7) along with the constraint a + \ > we get 

1 5-7V-2A; 
°< a+ 2 < N(N-l) - 



< 1. 



(3.7) 



(3.8) 



Let us consider the situation when k = 0. This corresponds to a class of eigenstates which 
includes the ground state of the system p], || and for N = 3 corresponds to the case of the zero 
angular momentum sector @. In this case we have 



< a 



1 5 — N 

2 < N(N-1 



(3.9) 
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Thus the allowed values of N which satisfy the above constraint are N = 3 and 4. These are 
the only cases where a negative energy bound state can appear in the system. For N = 3, the 
allowed ranges for a and g are given by— \ < a < — | and > g > — | respectively. Similarly 
for iV = 4, — |<a<— ^ and > g > — 1|. From the above discussion it follows that negative 
energy bound states may occur only when the 2-body coupling g in Eqn. (2.1) is attractive. 

Let us now consider the case when k > 0. From Eqn. (3.8) it follows that 5 — N — 2k > 0. 
There is no value of N > 2 which satisfies this constraint. Thus there are no negative energy 
bound states for k > 0. 

We therefore see that the original TV-body problem of Eqn. (2.2) admits negative energy 
bound states when k = and when iV = 3 or 4 with corresponding range of the coupling g as 
described above. Note that for k = 0, Pk{x) = 1. The full solution of Eqn. (2.2) in this case is 
given by 

tJj = By a+ ^r~^ +u \yfEr)^(J v (i^Er) - e l7TU J^ v {i^Er)). (3.10) 
The corresponding energy has the value £ = —E where E is given by Eqn. (3.6). 
We end this Section with the following observations : 

1) In the discussion above we have seen that the constraint a + | > puts strong restrictions 
on the allowed values of N. This constraint is needed to keep the wavefunction in Eqn. (2.3) 
well behaved when Xi —>■ Xj. It is possible to consider a situation where there is a natural cutoff 
in the system which prevents any two particles coordinates from coinciding with each other. In 
the presence of a natural cutoff, the wavefunction in Eqn. (2.3) would be well behaved even if 
a + \ < 0. In such cases it seems likely that negative energy bound states would exist for all 
values of N. 

2) When g = — |, or equivalently v = 0, Eqn. (2.12) admits an infinite number of bound states 
f7|, §. In this case a satisfies the condition a+ | = f^zfy > which together with the constraint 
a + | > precludes any value of iV > 2. However, when N = 3 and k = 0, the above condition 
yields a = — | which corresponds to a system of free particles with g = 0. Eliminating the 
centre-of-mass degree of freedom when N = 3 and g = 0, and using Jacobi coordinates yi, y 2 
|l|, Eqn. (2.2) can be expressed as 

/a 2 d 2 \ 

This represents the eigenvalue problem of two free particles in one dimension, or equivalently, 
one free particle in two dimension. Using polar coordinates and restricting to the zero angular 
momentum sector, the radial part of Eqn. (3.11) can be transformed in the form of Eqn. (2.12) 
with g = — |. Following 0, the solutions of Eqn. (3. 11) can thus be obtained as 

i> n = K (^£~ n r), £ n = -e^ l -^ cot ^ (3.12) 

where r = \Jy\ + y%, K (y/—£ n r) is the modified Bessel function and z represents the self- 
adjoint extension parameter. Thus we find that the zero angular momentum sector of this 
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system in Eqn. (3.11) admits an infinite number of negative energy bound states. This may 
seem surprising as the operator in Eqn. (3.11) is normally assumed to be positive. This is 
however not true in the presence of self-adjoint extensions which may break the positivity [ |15[] . 

3)Finally, when g < — |, the eigenvalue problem of Eqn. (2.2) does not have a physically 
meaningful solution as the the ground state energy in this case is unbounded from below [ |i~3|j . 
The effective Hamiltonian H could still be made self-adjoint through appropriate choice of a 
domain. This is however not sufficient for the ground state energy to be well defined. The 
2-body problem in this region of the coupling has been analyzed using renormalization group 
techniques ||16|| . Investigation of the iV-body problem in this case is still open. 



4. Connection to Virasoro Algebra 



The Hamiltonian of a 2-body system with inverse square potential is known to be related 
to conformal symmetry and Virasoro algebra [p] , [P7|| . The form of H in Eqn. (2.13) suggests 
that the effective Hamiltonian for iV-body system should also have similar property. Consider 
the Virasoro generators L n = —x n+1 -£; and shift operators P m = where ro, m E Z. In terms 
of L_i and Pi, H can be factorized as 

H = (-L-! + + nPx) (4.1) 

where \i = \ + v |ll| . Commutators of H with L_i and Pi leads to other elements of the algebra 



and continuing this process produces all the generators JTTJ. In particular, the commutator of 
H with a general element L m gives 

[L m , H] = 2/i(/i - l)P 2 - m ~ (m + \){L_ x L m _i + L m _iL_i). (4.2) 

The right hand side of Eqn. (4.2) contains products of the Virasoro generators. While these 
product operators are not elements of the Virasoro algebra, they do belong to the corresponding 
enveloping algebra. Thus H is contained in the enveloping algebra of the Virasoro algebra and 



has a well defined action on the associated representation space flT2" |. 

The above discussion is valid for arbitrary values of g. We now focus our attention on 
the bound states described by Eqn. (3.3) which are obtained when — | < g < |. If the 
effective Hamiltonian H is indeed associated to a Virasoro algebra, it might be expected that 
the wavefunctions exhibit a suitable scaling behaviour. The wavefunction x( r ) i n Eqn. (3.3) 
however has no clear scaling properties. This is due to the fact that in the presence of the 
self-adjoint extension, the scaling operator in general does not leave the domain of H invariant 



10| . Consider now the behaviour of the wavefunction when all the particles in the system are 



close to each other. In this limit of r — > 0, the behaviour of the wavefunction x( r ) is given by 
Eqn. (3.5). From this expression we see that for small values of r, the leading scaling behaviour 
of the x( r ) is governed by r~ u+ ^ and the wavefunction exhibits approximate scaling behaviour 
in this limit. 
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5. Conclusion 



In this paper we have investigated the existence of bound states in a one-dimensional N- 
body quantum system with inverse square interaction. We have shown that this system admits 
negative energy bound states when either 1) N = 3, k = and > g > — | or 2) iV = 4, k = 
and > g > — Moreover, there exists only a single such state for each value of allowed 
value of g. To our knowledge, this is the first demonstration of the existence of negative energy 
bound states in this system. 

These bound states have been obtained by analyzing the effective Hamiltonian H. The 
existence of these negative energy states is related to the self-adjoint extensions of H. Each 
choices of the self-adjoint parameter z leads to a different domain on which H acts. Since the 
boundary condition associated with a given physical situation is encoded in the choice of the 
domain, the self-adjoint parameter z is directly related to the physics of the problem. This is 
clearly seen by the fact that different values of z lead to inequivalent expressions for the bound 
state energy (cf. Eqn. (3.6)). We therefore conclude that the negative energy bound states are 
found only when most general boundary conditions are considered. 

The effective Hamiltonian H is contained in the enveloping algebra of the Virasoro algebra 
and has a well defined action on the associated representation space. This may seem to imply the 
existence of scaling symmetry in the problem which should be reflected in the properties of the 
bound states. The self-adjoint extensions of H in general break the scaling symmetry. However, 
when all the particles in the system are close to each other, the bound state wavefunction in 
Eqn. (3.5) exhibits an approximate scaling behaviour. This appears to be the only remnant of 
the scaling symmetry at the quantum level. 

It would be interesting to consider the system described here in the presence of a natural 
cutoff which may lead to bound states for all values of N. 
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